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Abstract 

A time dependent variational approach is used to derive the equations of motion for the Xcj)^ model. 
The simultaneous evolution of the quantum fluctuations and of the classical part of the field is considered 
in a lattice of 1+1 dimensions. Different initial conditions corresponding to non equilibrium situations 
are considered and evolved in time. Such high energy localized configurations expand in the lattice by 
"bumps" which may change with time. The quantum fluctuations make the peaks be smoother and the 
expansion faster. 

PACS numbers: 02.60.Nm; 03.65.Db; 07.70.+k; 05.50.+q; 11.30.Qc; ll.15.Tk; 11.90.+t. 
Key-words: quantum fluctuations, condensate, hydro-dynamical expansion, non-equilibrium initial con- 
ditions. 

IFUSP- /2000. 



*email:braghin@if. usp.br 
^email:navarra@if. usp.br 







1 Introduction 



In relativistic heavy ion collisions, such as those performed now at RHIC and CERN, very dense and nearly 
baryon free systems are formed in the central rapidity region. These excited systems decay and produce 
a large number of secondary particles. In view of the large densities and multiplicities, with or without 
quark gluon plasma formation, it is very likely that these intermediate systems will exhibit some collective 
behavior. It is therefore reasonable to treat their expansion with hydrodynamics ||l|, ^, |3| and it is interesting 
to investigate the connection between the hydro-dynamical and the field theoretical approaches, see for 
example Q. In fact, the fluid approach has been taken so seriously that irregularities in the hydrodynamic 
flow pattern have been proposed as possible signals of a deconfined phase ||5|. 

Given the physical appeal of hydrodynamics, many of its aspects and underlying assumptions have been 
investigated during the last ten years. Some of the studied topics are thermalization, initial conditions 
for the fluid expansion, the equation of state (and possible phase transitions), the freeze out mechanism 
and final state interactions. The direct correspondence of the hydro-dynamical description with the field 
theoretical one is thus highly appropriate. 

In this context, at least two questions arise which have not been addressed so far and are discussed in 
the present article. The first is: i) which is the role played by the quantum field nature of this dense matter 
in the expansion? In other words, which mistake are we doing when we use classical hydrodynamics and 
neglect quantum effects? The second question is: ii) if the expanding matter is merged in a condensate, 
which is the effect of this condensate on the expansion? 

Whereas solving these equations is technically very complicated in full hydrodynamics, the interplay 
between classical background and quantum fluctuations can be studied in simpler scenarios. In this work 
we discuss the space-time evolution of the A(/>^ model in 1 -|- 1 dimensions at zero temperature. Solving non- 
perturbatively and self-consistently the equations of motion of the theory, we can follow the expansion of a 
highly energetic system composed by a self interacting scalar field. In particular, we can separately investi- 
gate the classical and quantum components of the system and estimate the role played by the condensate 
during the expansion. 

In spite of the enormous differences between this simple model and the realistic descriptions of the 
"fireballs" formed at RHIC, we believe that our study can be of interest and give us some relevant insights 
into the real physical situation. 

Due to the extreme complexity of realistic theories, such as QCD, one usually is obliged to consider 
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effective models wliich respect the major properties of tlie fundamental theory. The Xcj)"^ model shares with 
QCD the properties of asymptotic freedom and spontaneous symmetry breaking (SSB) in the frame of the 
Gaussian approximation. Due to the SSB, the field is decomposed into a classical and a quantum part. At 
the same time it may also represent the mesonic sector of the linear sigma model, which is often used to 
describe the dynamics of a gas of pions. In the context of cosmology the scalar field of the model may be 
considered as the relevant degree of freedom for inflationary models Q. 

Some progress has been achieved in the last decade in understanding the dynamics of the A(/)^ model. 
In particular, lattice calculations of the time evolution of out of equilibrium configurations, both at zero 
and finite temperatures, in 1+1 and 3+1 dimensions in continuous and discretized space-times have been 
performed 0, |, g H, 0, [ll . 



In this work the will be investigated in the Gaussian approach with the formalism discussed in Refs. 
. It is suitable for non homogeneous configurations of the fluctuations (the quantum part of the field) as 
well as of the condensate (as the classical part of the field (j){t) will be called). They can be thought of as 
two interacting "liquids" with continuous energy transfer and this view is under investigation. The text is 
organized as follows. In section II the Gaussian formalism for out of equilibrium systems is briefly outlined 
and the equations of movement are derived. In section III we discuss the numerical method employed here 
and developed in . The numerical results are presented in section IV for several different out of equilibrium 
situations and values of the coupling constant. The results and some perspectives are summarized in the 
final section. 

2 Time dependent Gaussian approximation to model 

The Hamiltonian density for a scalar field 4> with bare mass ttzq and coupling constant A is: 

H = \ (vr2(x) + {Vmf + ml<p\x) + ^<^'(x)) , (1) 

where the action of operators (p and vr in functional Schroedinger representation over a wave functional 
$ [(f){x)] =< (/'(x)l^ [(f)] > is given by: 

>=m\^[<pm >, 

7r|^' [(/.(f)] >= -i6/6(l){x)\'^ [(/)(f)] > . 
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For our variational calculations analyzed below, in the Schroedinger picture the wave functional evolves like 
the Schroedinger equation 

i^^^ [0(f)] = im] ■ (3) 

This is, thus, a non covariant formalism suitable for time dependent problems. 
In the Gaussian approximation at zero temperature ^ is parametrized by: 

^ W^)] = A^exp J dxdy6(j){x) (G'^{x, y) + iY.{x, y)) 5(/>(y) + i j dfyf , (4) 

Where 6(p{x,t) = 4>{x) — (f){x,t); the normalization is N, the variational parameters are the condensate 
4>{x,t) =< ^'|(^|^' > and its conjugated variable Tt{x,t) =< ^'|7r|^ >; quantum fluctuations are represented 
by the width of the Gaussian G{x,y,t) =< ^'|0(x)0(y)|^' > and its conjugate variable T,{x,y,t). 

In variational time dependent calculations we have to choose an action to be minimized in order to 



obtain the equations of motion. We take the well known Dirac action from [13|: 

d 



1 = j dt<^\ -H^\^ > . (5) 
In order to calculate it, we take the mean value of an operator O given by: 

< *|d|^ >= /p[0]^'*d* (6) 



By means of this average procedure the energy density p =< > can be calculated. Its temporal 

evolution will be investigated in section 4. 

Variations with respect to the variational parameters and their conjugated yield the following equations 
of motion (any repeated spatial index means integration over that variable): 

dtG{x, y) = 2 (G(x, z)^{z, y) + S(x, z)G{z, y)) , 



5E(_x,y) 

dt4>ix) = -7f(x), 



dMx,y) = (^2S(x,z)S(z,y) - -G-\x,z)G-'{z,y)j + (^^^ + -<A(x) ^ 



6Tr{x) 



SI b 

dtTt{x) = r{x,y)(t){y) + -(p{x) 





Where r(x, y) = — A + [ttiq + -2G{x, x)j 6{x — y). In this approximation the interaction term A(/> becomes 
quadratic, i.e., it contributes to a self consistent mass. These equations were generalized for the out of 
thermo-dynamical equilibrium using different methods in 0, 0. 
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Using the Gaussian ansatz (P) in the symmetric phase {(j) = 0) we need only two initial conditions for 
the temporal evolution of these equations, G{t = 0) and G{t = 0), which is proportional to its imaginary 
part S. An inspection of the equations of motion shows that for (j){t = 0) = (j){t = 0) = the variables (j) 
and vf will not change, i.e., the classical part of the field will remain constant (zero) during all times in the 
frame of this approach. 

Several numerical works have been done for the above equations of movement. The choice of initial 
conditions is entirely subordinate to the Gaussian approximation. If it were not Gaussian, we would have 
to consider three conditions instead of two [^] . The analysis of these equations by Boyanovsky et al and by 

shows that initial conditions (for homogeneous G and (f) ) are crucial for the time interval in which the 
system evolves towards the minimum and for the speed of the field evolution. 



3 Numerical method 

The numerical method for the temporal evolution of the initial conditions corresponding to the system 
described above used in this paper was developed in j^, This section is a brief review. One defines a 
generalized density matrix in a lattice as: 

/ 



(8) 



pi J 

where the density matrices (mean values) are given by: pi,j = \ < CLiXij + a- aj > which is hermitian and 
Kij = — < aittj > is symmetric, using the creation and annihilation operators. These operators can be 
written in a lattice with mesh size Ax in d spatial dimensions as: 

a{j) = -^{(/>(i)(A:E)^ +m(j)(Ax)^} (9) 
= -^{</)(j)(Ax)^ -i7r(j)(Ax)'^^} (10) 

Some of the calculated mean values in terms of the matrix elements of the above matrix are given by: 

G{i,j) =< (l){i)(t>{j) >= \^_^ ^e{p{i,j) - K{i,j)) 

[llx)^ (11) 

F(i,j) =< UmU) >= j^^Mpihj)+<iJ)) 
where F{i,j) = G{i,j)-^ /A + 4^{i,k)G{k,l)^{l,i). 
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The temporal evolution is governed by the Hartree (Fock) Bogoliubov energy, which can be parametrized 
in the following form: 

/ 



Iff 

2 ''^ ~ 6R 



\ -Dij -Wij 



(12) 



where the above matrices are given in terms of the parameters of the model 
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Wi,j = + K + >^Gi,i/'^ + i^xy (13) 

Dij = ^(-Aij + {ml + XGi,i/2 - (Ax)-^)^. .). (14) 
With these matrices, one can check that the Liouville-von Neumann condition is satisfied: 

iRij = [Hik,Rkj] ■ (15) 

These equations are equivalent to those obtained from the Gaussian approximation of the previous section 
[^. The time-evolution of the generalized density matrix can be then performed for given initial conditions. 

4 Numerical results 

In this section we will report results of numerical calculations for the equations of motion studied above in a 
lattice of 100 points with spacing Ax = 0.1 fm. By varying Ax ^ it was possible to assure the reliability 
of the results (dynamics does not change in this limit). The lattice spacing is always to be much smaller 
than the correlation length (^ = 1/mphys) which is quite smaller than the lattice size (10 fm). The first 
initial condition for the classical part of the field is given by: 

^{t = 0)=4>o tanh^ iy^) ' <^ = ^ = 0, (16) 

where xq is the center of the lattice and the coordinate x is discretized. This configuration is not stable 
since in a small region of the lattice the field has values different from the vacuum. 

The field configuration corresponding to ( p^ ) is shown in Figure 1 with a solid line for the model with a 
strong coupling constant A = 600/m~^ and physical mass /i = IQOMeV . In terms of energy it corresponds to 
a bubble of high energy density in the vacuum which can be interpreted as an in medium effect. This initial 
condition is plugged into the system of equations (0) and numerically evolved with the method described in 
the last section. Knowing the field configuration at all times we can compute the energy density distribution 
at all time steps for different scenarios. Periodic boundary conditions are considered in this work. 
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In Figure 2a we show the energy density distribution of the classical part of the field alone (which we 
call p) along the x direction for different times t = Q (thick solid), =0.1 (solid), = 1 (thin solid), 2 (thick 
dotted) and 3.5 (dotted) fm. As it can be seen, the energy density excess with relation to the vacuum value 
is distributed among the lattice by means of "waves". The expansion of the initially high energy density 
proceeds with two peaks for each side of the lattice. This feature may be expected for the massless Klein- 
Gordon equation which has a sinoidal solution due to the gradient terms. The energy density is proportional 
to the square of the field and consequently we obtain this "two bumps" structure. 

In Figure 2b the same initial conditions for the condensate are evolved in time with the contribution 
from the quantum fluctuations, i.e., considering all the four equations (^. For the quantum part along this 
work, it was assumed that, at the initial time, they are at the vacuum value, i.e.: 

G(r, t = 0) = Go S oc G(r, t = 0) = 0, (17) 

where Go is the value of the quantum fluctuations in the vacuum of the asymmetric phase (0 ^ 0) obtained 
from a gap equation. In this case the gap equation is given by: 

/x2 = m2 + ^Go(^2) + ^<Ag (18) 

The value of Go can be fixed by the physical mass p. There is, now, energy transfer from the classical part 
to the quantum fluctuations and back. Comparing figures 2a and 2b we can clearly see that the inclusion 
of quantum fluctuations accelerates the expansion. This is especially visible a,t t = 3.5 fm and even more at 
higher times. The effect is small but unambiguous. We observe an inversion in the heights of the two bumps 
that characterize the energy distribution, the outer bump becoming more pronounced in the quantum case. 
This is an indication that, quantum fluctuations increase the pressure in spite of the isoentropic character 
of the approximation. Furthermore, the energy density configuration is clearly smoother than in the purely 
classical case. It is worth to remember once more that the present approach is not yet completely "self- 
consistent" in the sense that for given initial "non-equilibrium" conditions the model should reflect this fact 
including, for instance, the coupling constants and quantum fluctuations. This is under study []T4| . 

In order to investigate the sensitivity of the conclusion found in the last paragraph to the shape of the 
initial conditions, we repeat the same calculations starting from the following field profile: 

^(i = 0) = <Ao {l + 0.3 sech^ (^^) } • (19) 

The corresponding field configuration is shown in Figure 1 with a dashed line. In this case we use a different 
coupling constant A = 60 fm~^. This choice changes the global normalization of the field and energy, 
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but it does not change the dynamics considerably. As a matter of fact, in the classical level, there is a 
scale invariance in the dynamics with relation to changes in such variable. Taking into account quantum 
fluctuations this scale invariance is broken (slightly for the range of parameters considered in this work). 
Evolving only the classical part of the field in the same way as before we obtain Figure 3a, which is very 
similar to Fig. 2a. Switching on the quantum fluctuations leads to results of Fig. 3b. Comparing Figs. 
3a and 3b we observe the same effect already seen in Fig. 2, i.e., the inclusion of quantum fluctuations 
accelerates the expansion. 

In the asymmetric phase the effective potential has a Mexican hat like form and in all cases presented 
in this work the total energy is smaller than the "barrier" (located at the origin of the (f) space) which 
separates the two minima of the potential. In the cases shown in Figures 2b and 3b there is a tunneling 
of the condensate from one side of the potential to the other. For some kinds of initial condition this was 



already shown in This effect disappears if the quantum fluctuations are switched off but remain when 
the lattice spacing is reduced as discussed below. 

The presence of quantum fluctuations make the energy density to be amplified in the earlier times and 
to expand faster. This amplification is stronger in the case shown in figure 3b. As time goes the energy 
density "waves" tend to decrease. Moreover, we notice small regions with energy density lower than the 
true vacuum of the corresponding phase in the beginning of the expansion, just after the energy density 
amplification happens. 

The total energy of the system is conserved during the evolution. We want to stress that the inclusion of 
quantum fluctuations changes the overall energy normalization of the system. Therefore, the initial energy 
profile is similar for all figures but the normalization is different as can be seen in the pictures. 

The equations of movement (0) exhibit Ultra Violet (UV) divergences in a continuum space-time which 
are exactly the same as those present in the GAP equation for the vacuum of the model. They require mass 
renormalization (and also coupling constant renormalization depending on the space-time dimension) which 



is exactly the same for the static and time dependent cases. In 11] analytical solutions for these equations 
of movement were found for a particular kind of initial conditions. They exhibit the UV divergences only 
at the initial time (t=0). As this divergence is eliminated by the renormalization (of the mass and/or 
coupling) the temporal evolution remains unchanged. This is also true in the lattice: the deviation of the 
system with relation to the vacuum state determines the temporal evolution. For a variation of the lattice 
spacing, i.e., as Ax (for low values as 0.02 fm were checked) the dynamics does not change considerably 
(with the corresponding mass rescaling). In fact, the continuum limit must be the same for the static and 
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time-dependent cases. The deviations from this (renormahzed) state are always finite and are not altered 
by vacuum redefinition (renormalization) . Mass renormalization changes the normalization of the energy 
density being, therefore, not important for the dynamical evolution. 

In fact, the difference between the classical system and the classical plus quantum system is, in part, 
related to a scale invariance. The classical level exhibits scale invariance but this is not present at the 



quantum level. These effects will be explicitly addressed in detail elsewhere [14|. 

Our results point out to the following simple picture: without quantum fluctuations a given localized 
spatial configuration of the field (j) expands, simply converting "potential" into "kinetic" energy. The inclu- 
sion of these fluctuations accelerates the expansion. This effect may be caused by the interaction between 
the particles and the condensate, by the tunneling between the two vacua or by these both aspects. The 
questions addressed in this work are related to those discussed in Ref. |T[. In that work, the expansion was 
driven by the effective potential, (which contains a A</>^ part) plus the contribution of a gas of pions. The 
former and latter contributions correspond respectively to our classical and quantum contributions. The 
total pressure is the sum of both the effective potential and the pion gas component and therefore also in 
that model the inclusion of quantum fluctuations increases the total pressure and the expansion rate. These 
conclusions also are in agreement with studies for the effects of quantum fluctuations in the Inflationary 
scenario where there is "acceleration" of the dynamics, see for example |1C, 17|. 



5 Summary 

We have analyzed the temporal evolution of non-homogeneous configurations of the model considering 
two different approaches: the classical equations of motion and the equations of motion in the framework 
of the Gaussian approach for the quantum fluctuations in a one-dimensional lattice. We have been able to 
study the influence of the quantum fluctuations on the classical field dynamics and vice versa. Quantum 
fluctuations make the energy distribution smoother. Besides that, we conclude that the quantum fluctuations 
accelerate the expansion of a non homogeneous configuration of the classical part of the field after having 
amplified the energy density. We noted a tunneling of the condensate from one side of the "Mexican-hat" 
potential to the other when quantum fluctuations are considered depending on the initial conditions. 

We would like to emphasize that in the limit of vanishing lattice spacing the dynamics is not qualitatively 
modified: the renormalization acts as to change the field and therefore energy density normalization but 
not the dynamics. We have been concerned mainly with short interval evolution, but the evolution of 
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thermal degrees of freedom will be analyzed elsewhere. Finally with the inclusion of temperature, we hope 
to establish a closer connection between our system and a hydro-dynamically expanding fireball. 
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Figure Captions 



Fig. 1 Initial field configuration. Solid line corresponds to (|T^) and the dashed line line to (p^). 

Fig. 2 a) Evolution of the energy distribution of the condensate at different times; b) the same as a) with 
the inclusion of quantum fluctuations. The initial field configuration is given by (|l6|). 

Fig. 3 The same as Fig. 2 for the initial configuration (19). 
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Figure 1 
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Figure 3a 
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